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ABSTRACT 

We define a scalar measure of the local expansion rate based on how 
astronomers determine the Hubble constant. Our observable is the inverse 
conformal dAlembertian acting on a unit "standard candle." Because this 
quantity is an integral over the past lightcone of the observation point it pro- 
vides a manifestly causal and covariant technique for averaging over small 
fluctuations. For an exactly homogeneous and isotropic spacetime our scalar 
gives minus one half times the inverse square of the Hubble parameter. Our 
proposal is that it be assigned this meaning generally and that it be em- 
ployed to decide the issue of whether or not there is a significant quantum 
gravitational back-reaction on inflation. Several techniques are discussed for 
promoting the scalar to a full invariant by giving a geometrical description 
for the point of observation. We work out an explicit formalism for evaluat- 
ing the invariant in perturbation theory. The results for two simple models 
are presented in subsequent papers. 
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1 Introduction 



Quantum gravitational back-reaction offers an attractive model of cosmology. 
The idea is that there is no fine tuning of the cosmological constant, A, 
or of scalar potentials. In fact there need not be any scalars. Inflation 
begins in the early universe because A is positive and not unnaturally small. 
Inflation eventually ends due to the accumulation of gravitational attraction 
between long wavelength virtual gravitons which are pulled apart by the rapid 
expansion of spacetime. Inflation persists for many e-foldings because gravity 
is a weak interaction, even at typical inflationary scales, and it requires an 
enormous accumulation of gravitational potential to overcome this. Since the 
process is infrared it can be studied reliably using quantum general relativity, 
without regard to the ultraviolet problem [§]. Because the model has only a 
single free parameter — GA, where G is Newton's constant — it can be used 
to make unique and testable predictions |§. 

The physical mechanism of back-reaction requires quanta which are mass- 
less on the scale of inflation but not classically conformally invariant. This 
rules out competition from most ordinary matter, but it does allow an effect 
from light, minimally coupled scalars. It has been suggested that significant 
back-reaction can occur in scalar-driven inflation, even at one loop 0, 0]. 
It has also been proposed that scalar self-interactions can give a significant 
back-reaction at higher loops in A-driven inflation ||. All these models in- 
volve fine tuning to keep the scalar light compared with the scale of inflation, 
so they are probably not relevant to phenomenology. However, scalars have 
the great advantage of being comparatively simpler to study than gravitons. 

With any model of back-reaction one encounters the problem of reliably 
inferring its impact on the cosmological expansion rate. For a perfectly ho- 
mogeneous and isotropic geometry one would compute the expansion rate by 
transforming to co-moving coordinates, reading off the scale factor, and then 
taking its logarithmic time derivative. But back-reaction derives from the 
gravitational response to quantum fluctuations, and these break homogene- 
ity and isotropy. The notion of a cosmological expansion rate must obviously 
have a reasonable generalization since the current universe is not perfectly 
homogeneous and isotropic, yet astronomers mean something by measuring 
the Hubble constant. However, it is not so clear how to represent this ob- 
servable in terms of quantum gravitational operators. 

Previous studies of back-reaction have tried to resolve this problem by 
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averaging over fluctuations to produce an effective geometry which is homo- 
geneous and isotropic. Then the cosmological expansion rate is computed 
from this effective geometry in the usual way. In one method the averaging 
is accomplished by taking the expectation value of the gauge fixed metric in 
the presence of a state which is homogeneous and isotropic [§, ||, |7|, || . Then 
the expectation value of the metric must be homogeneous and isotropic even 
though it is the average over quantum fluctuations which are not. The other 
technique is to enforce homogeneity and isotropy by spatially averaging the 
gauge fixed metric over a surface of simultaneity [[|, ||. 

Serious objections have been raised to both techniques. Unruh dislikes 
using the gauge-fixed metric |§, either in an expectation value or in a spa- 
tial average. He argues that certain variations of the gauge fixing condi- 
tion change the expectation value (or spatial average) of the metric in ways 
which cannot be subsumed into a coordinate transformation. Unruh there- 
fore maintains that even forming the expectation value (or spatial average) 
of the metric into coordinate invariant quantities does not purge these quan- 
tities of gauge dependence. He would prefer that back-reaction be studied 
with an operator which is itself an invariant, before taking the expectation 
value. He also disbelieves that averaging over a surface of simultaneity can 
be relevant to what a local observer perceives. 

A different objection has been raised by Linde. He is willing to use the 
gauge fixed metric — and both men accept the validity of quantum field 
theory in determining the time evolution of the Heisenberg field operators. 
However, Linde suspects that inferring back-reaction with expectation values 
invites a Schrodinger Cat paradox. This is because inflationary particle pro- 
duction leaves the long wavelength modes in highly squeezed states whose 
behavior is essentially classical. No matter what Heisenberg operator is used 
to measure the cosmological expansion rate, Linde would prefer to stochasti- 



cally [|10] sample its probability distribution rather than take its expectation 
value. 

The present work is an attempt to address the preceding objections. To 
avoid potential problems from using the gauge fixed metric we propose to 
infer the local expansion rate instead from the functional inverse of the con- 
formal d'Alembertian: 



=d» (V^gg^) - 1 -r . (i) 
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This operator, acting on a unit "standard candle" , 



A[g](x) = ^-1, (2) 

averages over the past lightcone, as astronomers do when compiling a Hub- 
ble diagram. In the slow roll approximation the observable gives —\H~ 2 for 
an arbitrary homogeneous and isotropic universe. It is therefore a reason- 
able candidate for measuring the local expansion rate when the universe is 
not precisely homogeneous and isotropic. And it is a scalar function of the 
observation point x M . 

Nothing can be done about the noninvariance associated with the fixed 
initial value surface upon which the Heisenberg state is defined. However, 
invariance under the subclass of transformations which preserve the initial 
value surface can be achieved by geometrically specifying the point at which 
A[g] is observed. In scalar-driven inflation this can be done by defining 
zero-shift surfaces of simultaneity so that the quantum inflaton agrees with 
its classical value. (Using these coordinates was Unruh's suggestion.) In 
more general models one can build invariant surfaces of simultaneity using 
the inverse minimally coupled d'Alembertian. The expectation value of the 
resulting invariant can then be evaluated, or else its probability distribution 
can be sampled stochastically. 

In Section 2 we motivate the scalar and show that it has the proper cor- 
respondence limit for exactly homogeneous and isotropic geometries. Section 
3 discusses the corrections needed to geometrically specify the observation 
point. In Section 4 we expand the scalar in powers of the metric fluctua- 
tions. Section 5 concerns the retarded Green's functions which appear in 
this expansion. We discuss a somewhat more complicated but considerably 
sharper observable in Section 6. Our conclusions comprise Section 7. Two 
subsequent papers give the results of applying the observable to models of 
scalar-driven [O] and A-driven [12] inflation. 



2 Motivating the scalar 

Since we are interested in the effect of back-reaction on inflation it is rea- 
sonable to consider perturbations about a background geometry which is 
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homogeneous, isotropic and spatially flat, 

ds 2 = -dt 2 + e 2h{t) dx ■ dx = a 2 (r]) {-drf + dx ■ dx) . (3) 

There is general agreement that the cosmological expansion rate for this 
background is H = b = a' /a 2 . Dots denote co-moving time derivatives while 
primes represent conformal time derivatives. We normalize the initial (t = 
or i] = rjj) scale factor to unity. 
The full metric has the form, 

9»AV, x) = a 2 (r])g^ u (r], x) = a 2 (r]) [rj^ + ki/j^tj, x)} , (4) 

where r]^ u is the spacelike Lorentz metric and k 2 = 16nG is the loop counting 
parameter of quantum gravity. Fluctuations reside in the pseudo-graviton 
field, ipnuir], x), whose indices are raised and lowered with the Lorentz metric. 
What we seek is a scalar functional of the metric which provides a reasonable 
extrapolation for how a localized observer would measure the cosmological 
expansion rate when ip^uirj^x) ^ 0. 

It is worth explaining why the Ricci scalar is not satisfactory. R(x) is 
certainly scalar, and it is closely related to the Hubble constant for the case 
of perfect homogeneity and isotropy, 

R — ► 12H 2 + Qa^H' = 6a~ 3 a" . (5) 

However, no local curvature invariant can account for the ability of observers 
to perceive the larger universe at cosmological distances by looking back 
along their past light cones. Einstein's equations set the Ricci scalar to 
— 8ttG times the trace of the stress tensor. This actually vanishes during a 
phase of radiation dominated expansion! Nor does the local value of R(x) 
have much to do with what an observer can see at cosmological distances. 
For example, even "empty" space within our solar system contains about 10 
Hydrogen atoms per cubic centimeter. Were we to infer the rate of cosmolog- 
ical expansion using R(x) the result would correspond to a Hubble constant 
about a hundred times larger than the actual value, 

8nGp ~ 3 • 10- 3 °i ~ ( 5 • 10 4 km ) 2 . (6) 

' s 2 v s-Mp C ; v 1 

We stress that there must be a reasonable solution to this problem because 
the current universe is not precisely homogeneous and isotropic, yet human 
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astronomers still claim to be able to measure the Hubble constant. It is 
instructive to review one of their simpler techniques. Consider light emitted 
at (771, and received at (770, x ). The observed quantities are the redshift 
z and the flux T . If the source luminosity C is known, these two quantities 
can be related under the assumption of perfect homogeneity and isotropy. 
Astronomers simply define the local Hubble constant so as to make the same 
relation true in the presence of fluctuations. Then they average over many 
sources. 

Deriving the relation between T and z is a standard exercise [|TJ|. Assum- 
ing perfect homogeneity and isotropy the physical distance between source 
and observer at time 770 would be, 

Ar = a \\x - xx || = 00(770 - Vi) ■ (7) 

The measured flux is the flat space formula, corrected for the redshifts of 
energy and rate, 



1 + zJ 4vrAr' 

One inverts this relation to solve for the product of 1 + z times Ar, which is 
known as the "luminosity distance," 



d^(l + ,)Ar = ^. (9) 

If both observer and source are at rest in conformal coordinates then the 
observed redshift would be, 

* = -rr-i- (10) 
0(771) 

Its relation to Ar comes from the scale factor's Taylor expansion, 

1 



1 - H Ar + ^(1 - q )H*Ar 2 + . . . 



0(771) = a 

where the current Hubble constant and deceleration parameter are, 



(11) 



H o = -2 , Qo = 1 jt ■ (I 2 ) 

"0 
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Inverting to solve for Ar gives, 



H Ar = z-^(l + q )z 2 + ... . (13) 

and multiplication by 1 + z results in the luminosity distance, 

H,d L = z+ l -{l-q )z 2 + ... . (14) 

Plotting z against for relatively small z gives a straight line whose inverse 
slope is the Hubble constant. 

It is not simple to identify invariants which represent the observed quan- 
tities z and T for an arbitrary metric. If one considers the transmission 
process in terms of individual photons then the redshift could be formulated 
as follows. Let us denote the worldlines of the emitter and observer as func- 
tions of their respective proper times by X£ m {r) and X^ hs (r). Recall that 
proper times are normalized to obey, 

g^(X(r))X a (r)X^r) = -l, (15) 

where dots stand for differentiation with respect to r. Now consider a photon 
which was emitted at proper time T\ and reaches the observer at proper time 
To. Of course the afhne parameter a of the photon's worldline Xp h (a) cannot 
be a proper time since the 4-velocity must be lightlike, 

g aP (X ph (a))X^ h (a)X^(a) = . (16) 

Given any function X^ h (a) which obeys ( |i~6"|) as it interpolates from X^ m (ri) = 
Xi to X^ bs (ro) = Xq, one makes a reparameterization of the afhne parameter 
(the new value of which we shall continue to call a) so as to enforce the 
geodesic equation, 

+ rv K h (a)x; h (v) = o . (17) 

The redshift experienced by such a photon is given by, 

1 + z = ffp^l)^ph(^l)^cm(>l) ^ lg j 

g^(x )X£ h ((To)X» bs (T ) 
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The flux is essentially the response, at the observer's location, to the 
(presumed known) source's current density J M (x). One begins by solving 
Maxwell's equations for the field strength tensor F^ u (x), 

Fpu" — J fx j Fa0;-y + F/3r,a + Fya;f3 = • (19) 

These equations are invariant under local conformal rescalings, 

Ffiv ¥ F'/j.u 1 Jfi ¥ ^ J (i > 9[iv * ^ 9[iv ; (20) 

and they can be solved in terms of something we shall call the conformal 
tensor d'Alembertian, 

^□f = ^ - 8^W V - b\R\ + 2J*7„ . (21) 

The solution is, 

Ffj, u — pu ^ — J {—Jp;a + Jcr-p) ■ (22) 

One gets the stress tensor from the field strength tensor, 

t,u = - \g»»g aP ) g pa F ap Fp a . (23) 

The Poynting vector is obtained by contracting the observer's 4- velocity into 
the electromagnetic stress tensor, 

Sp = T^(xo)X^(ro) . (24) 

And the measured flux is the norm of the Poynting vector in the orthogonal 
projection of the observer's metric, 

T 2 = S,S U {f + X&jZ*) ■ (25) 

Astronomers measure electromagnetic radiation because it is available to 
them, but this choice of observable complicates the metric dependence of 
the operators which represent their measurements. For example, the tensor 
character of Maxwell's field strength is why one has to invert the tensor 
conformal d'Alembertian in (22), rather than its simpler scalar cousin. It is 
also why the response field has to be squared in (E3p. Other complications 
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arise from the fact that the source luminosities are not precisely known, and 
that their distribution throughout space is not uniform. 

The preceding complications pose important limitations on observational 
astronomy but they need not restrict our choice of the operator with which to 
probe the theory. For us the really essential feature is to measure the response 
of some long range field to known sources distributed along the observer's 
past lightcone. We can retain this feature and vastly simplify our labor by 
observing a conformally coupled scalar, rather than a conformally coupled 
tensor. We can achieve a further simplification by taking the source to be a 
uniformly distributed monopole, rather than a sparse distribution of dipoles 
of varying strength. Then a single measurement of the scalar represents a full 
sky average and we can dispense with the complication of having to tabulate 
two quantities (z and J 7 ) for each source point. We call the scalar A and 
define it to obey the equation, 

O c A = 1 , (26) 

where d c is the conformal d'Alembertian ([!]). If we define the scalar and its 
first time derivative to vanish on the initial value surface the result is just 
the integral of the retarded conformal Green's function, 

A[g]{x) = 1-1. (27) 

It remains to show that *4,[g](:r) has the right correspondence limit for 
exact homogeneity and isotropy. In this case the conformal d'Alembertian 
reduces to the form, 

□ c — ► a~ 3 d 2 a , (28) 

where d 2 = rf^d^dy is the flat space d'Alembertian in conformal coordi- 
nates. The operator becomes even simpler acting on spatial constants. One 



consequence is that (26) can be solved by simple integration, 



Ao(rj,x) = -a- 1 ^) f V f ' drj"a 3 ( V ") = -e~ b ^ f dt'e^ f dt"e 2b ^ . 

Jr\i Jr\i JO Jo 

(29) 

It turns out that ( |29"D can be evaluated quite generally in what is known 



as the slow roll approximation. This is obeyed by all successful models of 
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inflation and it amounts to neglecting all higher co-moving time derivatives 
of the logarithmic scale factor b(t) with respect to the first, 



d N b 



dt N 



«(b) 



N 



V N > 2 



(30) 



Most operations of ordinary calculus can be done explicitly in the slow roll 
approximation. For example, the following trivial rearrangement, 



,26 




+ e 



26 



2b 2 



(31) 



allows us to express the initial integrand of (|29"D as a total derivative plus a 
term which is negligible in the slow roll approximation. It would be straight- 
forward to develop a series in slow roll corrections but the first is generally 
sufficient for the inflationary setting in which we wish to employ the new 
observable. With positive exponents and any significant amount of inflation 
it is also possible to ignore the lower limit, 



dt'e 2b ^ 



e 2b(t') 

2b(t') 



2b 2 [f] 



+ ... 



e 2b(t) 

2b(t) 



We can therefore apply the slow roll approximation to Eq. (|2 
the observable in terms of the Hubble constant, 



(32) 



to express 



Ao(r),x) 



1 



2b 2 {t) 



(33) 



In analogy to astronomical practice the local Hubble constant in the presence 
perturbations is defined so as to preserve this relation, 



A[g]{x) 



1 



2H 2 (x) 



(34) 



Although we are chiefly interested in applying the new observable during 
inflation it worth noting that the slow roll result (|34] ) is valid, up to a number 
of order one, for quite general geometries. For example, with general power 
law expansion the logarithmic scale factor and Hubble constant are, 



b{t) = s In 1 + 



b(t) = Hj(i + 



Hit 



(35) 
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where Hi is the initial Hubble constant and s is a constant. With this simple 
time dependence we can perform the integrals in fl29[) exactly, 



*power 



s 2 e -b(t) g 2 g-26(t) 

+ ( s _l)( s + 2 ) H] ~ ( s -l)( s + l)2H I b(t) ' (36) 

The second and third terms become insignificant at late times and the first 
rapidly approaches (|3^ ) for large s. Even for s = | the numerical factor is |. 

3 Fixing the observation point geometrically 

Even scalars depend upon the point at which they are observed. Part of 
this dependence is physical. The Heisenberg state is specified on a particular 
initial value surface and the geometrical relation of the observation point to 
this initial value surface can and should affect the result. The purpose of this 
section is to formulate the technology for imposing such a relation. 

Of course any method of describing points amounts to fixing a gauge 
however, there is an important distinction between ad hoc gauge conditions 
and those which exploit some special feature of the particular system under 
study. For example, if the system includes a Sun then it is geometrically 
meaningful to take this star's center as the spatial origin. It is therefore 
necessary to be as precise as possible about the nature of the system under 
study. 

The dynamical variables of our system include the metric g^x) and pos- 
sibly also a scalar inflaton field (f(x). Our goal is to compute the expectation 
value of the operator A[g] (or to stochastically sample its probability dis- 
tribution) in the presence of a Heisenberg state which we shall assume is 
homogeneous and isotropic. Since the only effective technique for making 
such a calculation is perturbation theory we shall also assume that the var- 
ious quantum field operators are perturbations on a background which is 
homogeneous and isotropic, 

(p(r},x) = (p (rj) +cj)(r),x) , (37) 
9ij,u(v,x) = a 2 (r])(r] flL , + nijj^fax)) . (38) 
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Because geometrically significant gauge conditions can involve nonlocal 
and nonlinear functionals of the fields, we wish to preserve the option of 
carrying out the calculation in more convenient gauge. A simple technique for 
accomplishing this is to define the observation point as the field-dependent 
coordinate transformation Y^tp, g](x) such that the transformed scalar (if 
there is one) and the transformed metric, 

<p'(x) = <p(Y(x)) , (39) 
gyp QY a 

CM - a^a^^(Y(x)) , (40) 

obey the geometrically significant gauge conditions. Then one can evaluate 
«4[<7](y(:c)) in any gauge and the result will be the same. 

The existence of a fixed initial value surface (at rj = rji) suggests that 
Y^ should be expressed as the composition of a temporal transformation 
i] — > t(i],x) followed by a purely spatial transformation x l — > x t { r )i 
Surfaces of simultaneity are defined by the condition t(i],x) = constant, 
while x t { f ]>x) traces out "the same" space point on the foliation of these 
surfaces. The full transformation would be, 

Y°(r],x) = t(t),x) , Y\ri,x) = (r(ri,x),x) . (41) 

The problem's homogeneity and isotropy implies that all space points are 
physically equivalent and we may as well use orthogonal projection to define 
"the same" space point. This amounts to the condition g' Qi = and hence, 

°-9-f9o k (v,x) + -^g-9 jk (v,x). (42) 
The relation can be simplified by multiplying with the inverse Jacobian, 

9 Mx) + ^9ij(v,X)=0- (43) 

Whereupon multiplication by the inverse 3-metric results in the following 
first order (but nonlinear) differential equation, 

<>v : -(g-'fgo^x), (44) 



dr] 

= -Ki/J i(v, X) + K 2 (fojipji) (v, X) - ^ {i>oki>kji>ji) (v, x) + ■ ■ ■ (45) 
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Making the obvious choice of initial condition gives an integral equation 
whose iteration to any order is straightforward, 

X l (v,x) = x*- f V da(g- i y g0j (a,x), (46) 

rv rv 
= x l — k I daip Qi (a,x) + k i da (ip jipji) (a,x) 

'vi J VI 

fv r a 
+k 2 I datp 0iJ (a,x) dpip 0j (p,x) + 0(/t 3 ) . (47) 

J VI J Vi 

Defining surfaces of simultaneity is less subtle for A-driven inflation than 
for its scalar-driven cousin. Without back-reaction the cosmological expan- 
sion rate is constant in A-driven inflation, so the effect is certainly real if one 
sees progressive slowing under any timelike foliation. In the scalar-driven 
case there is already slowing as the background scalar rolls down its poten- 
tial so one must be careful to compare the expansion rate with and without 
back-reaction at the same physical time. 

Since the value of the scalar determines the expansion rate without back- 
reaction it seems reasonable to define surfaces of simultaneity so that the full 
inflaton field agrees with its background value, 

<p(r(ri,x),x) = <p (rj) . (48) 

This can be solved perturbatively by first writing, 

r(r], x ) = 7] + 5r(rj, x ) , (49) 

and Taylor expanding, 

n=l U - n=0 n - 

Inverting results in an expansion for 5r in powers of the quantum scalar <fi 
and its derivatives (all evaluated at (r],x)), 

(f) 00' (p'n(f) 2 , . 

^ = --S + ^-Sr + ^)- ( 51 ) 

A-driven inflation can be included within the same scheme by employing 
a scalar functional of the metric with monotonic time dependence in place of 
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Figure 1: Invariant procedure to fix the observation point. 



the dynamical scalar. Perhaps the simplest of these "clock functions" makes 
use of the inverse minimally coupled d'Alembertian acting on the Ricci scalar, 



SV[g]{x) = -—R 



(52) 



We define surfaces of simultaneity so as to make the clock agree with its 
background value, just as relation (pE8|) does for scalar-driven inflation, 

jV(r(r],x),x)=jVo(v) • (53) 

Fig. 1 depicts the resulting foliation. 

To see that j\f\g)(x) is a good clock in perturbation theory note that 
jVq(i]) ~ \n(a(j])) in the slow roll approximation. This follows because the 
minimally coupled d'Alembertian takes the following form in a homogeneous, 
isotropic and spatially flat geometry, 



(54) 



When acting on functions of only time the inverse of □ reduces to, 

I rv rv' rt rt 

□ Jt]i Jf)i Jo 



1)1 



dt'e- 3b{t,) f dt"e 3b ^ . (55) 
Jo 
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With no perturbations we therefore have, 

ft rt' / Q 



lVo(v) = J dt'e' 3b ^ J dt"e 3b ^ (3b 2 (t") + ^6(i")) • (56) 



Making the slow roll approximation gives, 

■A/o(^) « f Q d1*e-w£df±(b(1?')eWn) , (57) 
« /* dt'fe(t') = &(t) . (58) 

JO 

Before closing the section we should comment that there is no problem, 
in perturbation theory, about evaluating an operator such as ^4[^](x) at a 
point such as Y^[(p, g](x) which is itself an operator. One merely expands in 
powers of the perturbatively small quantity Y^(x) — x^. Only a finite number 
of terms need be included to reach any fixed order in perturbation theory. 
Note that the various operator products should be time-ordered. This is 
because the functional integral of the ©-number functional A[g] (Y[ip, g](x)) 
is manifestly invariant as it is, and gives the expectation value of the time- 
ordered product of the corresponding operator. 

4 Pseudo-graviton expansion 

The purpose of this section is to expand ^4[^](x) in powers of the pseudo- 
graviton field ip^u(Vy £)■ This is most easily accomplished by first expressing 
□ c in terms of the conformally rescaled metric, 

9^u(v, x) = a -2 (77)5^(77, x) = rj^ + m/v(t?, x) . (59) 

We now write O c = a~ 3 Va, where V and its expansion are, 

1 „ / r— _.„.„ \ 1 



T> = ^pf* [yJ-9 9^9 V J -^R = d 2 + kPx + k 2 V 2 + ... . (60) 
The first two operators in the expansion are, 

v x = -r v d,o u + (~r a a + \^) (V v - t p P ) , (6i) 
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We remind the reader that pseudo-graviton indices are raised and lowered 
by the Lorentz metric, r]^. Other notational points are that the trace of 
the pseudo-graviton field is ip = ip p and that the second order, conformally 
rescaled Ricci scalar is, 



R 2 = r p (^V^- 2 ^7 

The next step is to factor d 2 out of T> 



- -A a ^^ a - + ^ Q/ V-« - - A t a t« ■ (63) 



V = d 2 (l + + ^ 2 P 2 + 0(k 3 )) . (64) 

Inverting D is now straightforward, 

111 11 1 11,1^.,. 
V = & ~ V^W 2 + ¥ KVl ^ KVl ¥ " 'd 2 'd 2 + 0{K ) ■ (65) 

All this implies the following expansion for the scalar observable, 

A[g] = a ~ 1 ^ fl3 = A + kAi + k 2 A 2 + 0(k 3 ) . (66) 
Ao was worked out at the end of Section 2. The next two terms are, 

Ax = -a-^V^a* , (67) 

^ = -a-'^V^ + a-'^V^V^a 3 . (68) 

It should be noted that the pseudo-graviton field is not free, nor are all 
of its components dynamical. The next step after this would be to expand 
ipfiuiVy m terms of the fundamental dynamical degrees of freedom, whatever 
they happen to be. This obviously depends upon selecting a particular model 



and must be postponed until this has been done |TT], [12 



5 Retarded Green's functions 

The pseudo-graviton expansion of the previous section results in a series of 
terms which involve the inverse differential operator 1/d 2 . The purpose of 
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this section is to precisely define the action of this operator. We also apply 
the slow roll approximation. 

The first task is easily accomplished. The retarded Green's function for 
the operator d 2 is well known, 



G(x; x') 



-5(Arj - Ax) 



(69) 



AttAx 

where Arj = rj — rf and Ax = ||5 — x! ||. Since the initial value surface is at 
77 — rji we define the result of acting 1/d 2 on an arbitrary function f(rj,x) 
as, 

rv . r „ ,5(Ar]-Ax) 



If 

d 2 \ 



(r),x) 



d V ' ld*x' ^ A -' f( V ',x> 



drj' Arj 



d 2 n 
Air 



f (r/, x + Ar]n) . 



(70) 
(71) 



When the function depends only upon time we can reach a form similar 
to that of Section 2. Making the substitution f(r],x) — ► Fiji) gives, 

1 



Q2 



(77, x) = — dr]'Ar]F(r] ) , 

j m 

rv rn' 
= - dr/' dri"F{ri") . 

JriT JriT 



(72) 
(73) 
(74) 



An important example is provided by the rightmost term for each of the *4 n 's 
- d~ 2 a 3 . We can explicitly evaluate these terms by making use of the slow 
roll approximation, 



d 2 



(77, x) 



2b 2 (t) 



(75) 



The fact that this depends only upon 77 allows one to simplify expressions in 
which derivatives act upon it. For example, the contribution from the first 
term in (^l|) is, 
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1 / 3 

-a '—Kyj 00 a 



d 2 



n 



a^(r]) I dT]'a 3 (r]')Arj 
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d 2 n 
Air 



ki/j 00 (77', x + Arin) . 



(76) 
(77) 



Further progress requires using the slow roll approximation. It turns 
out, as one evaluates the various factors of 1/d 2 from left to right, that 
the various integrands upon which they act are always dominated by the 
universal initial factor of a 3 . In typical gauges the pseudo-graviton field can 
grow at most like powers of ln(a). Although derivatives can sometimes result 
in a net loss of powers of the scale factor, they can never add such powers. 
Further, whenever even a single power of a is lost the contribution which 
finally results to *4.[g](x) is exponentially suppressed and hence irrelevant. It 
therefore suffices to consider terms of the form d~ 2 (a 3 f) for functions f(r), x) 
which grow less rapidly that a(n), 



d 2 



a 3 f 



jo J An 



(78) 



Note that A77 and 7/ are the following functions of t and t', 

A V = fdfe-WI , i = m+ f dfe-W'K (79) 
Jt> Jo 

Because An vanishes at t' — t a single partial integration fails to extract 
the leading order term in the slow roll approximation, 



C dt'e 2b(t,) Ar]f (n',x + Ann) 
Jo 
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Arjl 
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(80) 



The surface term vanishes at the upper limit and is exponentially suppressed 
at the lower limit. The remaining integrand is, 



dt' \ 2b ) 2b 2 2b 1 



Arjn ■ Vf 



51) 



The first term on the right is the original integrand times a term which is 
negligible in the slow roll approximation. Of the remaining terms only the 
one without the factor of A77 survives at the upper limit after another partial 
integration. Since each additional partial integration produces either a factor 
of b/b 2 or of e~ b the slow roll approximation of (|78|) is, 
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d 2 



a 3 f 



(r),x) 
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(82) 
(83) 
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The vast simplification inherent in (|83| ) derives from the fact that the 
response of a conformally coupled scalar at 77 = r] to that part of the source 
at 77 = rji redshifts as 0(770) / a (Vi) ■ Hence only the most recent sources 
contribute effectively. One can recapture the slow roll approximation much 
more simply by rewriting the differential equation fl26|) which defines *4.[g](x), 



Since the Ricci scalar is almost constant during inflation the leading order 
slow roll term — to all orders in the pseudo-graviton expansion — is con- 



6 A better observable 

Simple expressions are nice but the results of the previous section are too 
much of a good thing. To leading order in the slow roll approximation our 
new observable has turned out to be nothing more than —6 over the Ricci 
scalar! The next order terms disrupt this correspondence but still give ex- 
pressions which are local in the observation point. We criticized this sort of 
locality in section 2. A reasonable measure of the cosmo logical expansion 
rate should not be dominated by local fluctuations. Of course small fluctua- 
tions are, by definition, sub-dominant to the background, so one can still use 
*4.[g](x) to measure back-reaction during the first stages of inflation. How- 
ever, the defect of too much locality has a relatively simple fix which we shall 
present in this section. We shall also refine the observable so that it gives 
the Hubble constant exactly for any homogeneous and isotropic geometry, 
without recourse to the slow roll approximation. 

The effective locality of .4.[g](a;) derives from the fact that conformal 
scalars redshift like the inverse scale factor. During inflation the scale factor 
grows so rapidly that only the most recent sources matter much. A straight- 
forward way of avoiding this is by breaking conformal invariance. Suppose 
we measure a minimally coupled scalar B[g}(x), whose value and whose first 
derivative vanish on the initial value surface, and which is driven by a source 




(84) 



tained in the reduction, *4.[g](x) 



6/R(x). 



S(x), 



UB{x) = S(x) . 



(85) 
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For a homogeneous and isotropic geometry ([3]) the slow roll approximation 
results in an almost uniformly weighted average over comoving time, 

-S — - fdt'e 3 ^ fdt"e^S(t") (86) 
□ Jo Jo 



« -fdt'^l. (87) 

It remains to identify a suitable source. Note that for a homogeneous and 
isotropic geometry (^) the nonzero components of the Ricci tensor are, 

R 00 ^-3b-3b 2 , R l3 — > (b + 3b 2 ) 9ij . (88) 

The trace of the spatial part has the curious property of giving a total deriva- 
tive when multiplied by e 3b , 

3 (b(t) + 3b 2 (t)) e 3fc « = j (3b(t)e 2b ^) . (89) 

If the source S(x) reduces to minus one third times this spatial trace, for a 
homogeneous and isotropic geometry, then the minimally coupled scalar will 
reduce to the logarithmic scale factor exactly, 

B[g](x) — > jf* dt'e 3b{t,) £ dt"e 3b{t,,) (b(t") + 3b 2 (t")) = b(t) . (90) 

We can then obtain the Hubble constant by differentiation with respect to 
comoving time. 

The object we have just described is not quite a scalar because the spatial 
trace of the Ricci tensor is not. However, we can give the latter an invariant 
formulation by exploiting the technology of section 3 to define it in a special 
coordinate system which reduces to conformal coordinates for a homogeneous 
and isotropic geometry. With the transformation y M [g](x) we can define the 
spatial components of the metric and the Ricci tensor, 



9 '- {x) ^^^ 9 " AY{x)) ' ^ {x) = -^-dx7 R >" {Y{x)) ■ (91) 
The 3-curvature is just the inverse of the first contracted into the second 



K[g](x) = (g'-^R'Jx) , (92) 
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and the minimally coupled scalar is, 

B[g]{x) = ^(- l -n) . (93) 

Since conformal coordinates have zero shift we obviously wish to use the 
same spatial transformation (|46| ) as in section 3. The temporal transforma- 
tion requires a scalar clock function, the most uniformly applicable choice of 
which is V[g](a;), the invariant volume of the past lightcone as seen from the 
point x^ back to the inital value surface. Note that it must be a good clock 
generally, not just in perturbation theory, because the volume of the past 
lightcone increases monotonically under any timelike foliation. 

We define V[g](a;) as the invariant integral over all points which are con- 
nected to x' M by any future-directed, non-spacelike path.Q For a homogeneous 
and isotropic geometry this reduces to a single integral, 



V[g](x) — > V Q (v) = U I* drf&fflir, - r/) 3 . (94) 

6 Jri! 



As in section 3, we define surfaces of simultaneity to make this relation persist 
in the presence of perturbations, 

V(t( V ,x),x)=V (v). (95) 

We define the general conformal factor as the square root of the 00 com- 
ponent of the metric in these coordinates, 



Q[g](x) = JgUx), (96) 



dr 
dr] 



9oo(Y(x)) ~ [g 1 Y J 9oi9oj(y(x)) 



(97) 



Since the coordinates have zero shift, g' oi = 0, and it follows that Q is precisely 
the factor needed to scale from conformal to co-moving time. One possible 
definition for the Hubble constant is therefore the co-moving time derivative 
of the scalar B evaluated in these coordinates, 

H x \g\{x) = Sl- l [g]{x)^B[g] (Y[g](x)) . (98) 



1 The path need not be a geodesic, nor does it have to be the sole path which connects 
and x^. 
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Of course one might equally well base the observable on Q [g] (x) now that 
we have it, 

H 2 [g](x) = n-^Kx^ln (%](*)) . (99) 

One might instead employ the third root of the determinant of glAx). All of 
these are plausible measures for the cosmological expansion rate, all reduce 
exactly to the Hubble constant for homogeneous and isotropic geometries, 
and we anticipate that all will give the same result as regards the existence or 
non-existence of a significant back-reaction. We emphasize this multiplicity 
of plausible observables is as it should be because a similar situation exists 
in the many different methods by which astronomers attempt to measure 
Hubble constant. 



7 Discussion 

Reliably quantifying the effect of back-reaction on inflation poses a frustrat- 
ing paradox. The possibility of an effect derives from fluctuations in homo- 
geneity and isotropy, but these call into question precisely what is meant by 
the rate of cosmological expansion. Previous work has attempted to resolve 
the issue by averaging the gauge fixed metric, either over a surface of simul- 
taneity @, || or over the range of quantum fluctuations in a homogeneous 
and isotropic state 0, |(| |7|, || . It has been objected that neither technique is 
manifestly invariant, and also that the former procedure involves superposing 
data unavailable to a local observer on the surface of simultaneity ||. 

In section 2 we argued that both problems can be avoided by measuring 
the response of a noninteracting, conformally coupled scalar to a constant 
source. The scalar ^4.[g](:r) is a nonlocal functional of the metric which is ob- 
tained by superposing over the past light cone of just as astronomers do 
in measuring the Hubble constant. Its phenomenological interpretation also 
follows the standard practice in astronomy: we define the locally observed 



rate of cosmological expansion to bear the same relation — equation (|34j) - 
to *4.[g](:r) for a general metric as it does for a homogeneous and isotropic 
one. Of course the result will be a little different at different locations, just 
as we must expect the Hubble constant measured by human astronomers to 
disagree slightly with the value obtained from the different field of view avail- 
able to their opposite numbers in the Coma Cluster. But nearby observers 
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will tend to agree because their past lightcones largely overlap. 

It should be noted that we are not adding a conformal scalar to whatever 
model of inflation is being probed. The observable is only a functional of 
the metric used to pose invariant questions about the expansion rate; it does 
not change the dynamics of the model. Even if one insists that the scalar 
represents a sort of measuring device whose effect must be included, the 
strength of the constant source can still be adjusted so as to make this effect 
negligible. The required constant would simply appear on the right hand 
sides of both equations (|27|) and fl34[) , 



AWx) ^ 1-K ^ _0_ y {m) 

so that the magnitude of the scalar could be made arbitrarily small without 
affecting our determination of the local Hubble constant. 

One can either compute the expectation value of ^4[(7](x) — or else sto- 
chastically sample its probability distribution — in the presence of a Heisen- 
berg state which we assume to be homogeneous and isotropic. The passage 
from a scalar to an invariant can be achieved by geometrically fixing the ob- 
servation point relative to the initial value surface on which the Heisenberg 
state is defined. Because all points on the initial value surface are physically 
equivalent the problem reduces to orthogonally projecting between geomet- 
rically specified surfaces of simultaneity. Two definitions for such surfaces 
were presented in Section 3, along with perturbative expansions for the field 
dependent observation point Y^\(p, g](x) which can be used to fix the obser- 
vation point when working an arbitrary gauge. 

Sections 4 and 5 developed the general machinery necessary to evaluate 
the new observable perturbatively. We emphasize that these computations 
are imminently doable in the slow roll approximation. It remains to apply the 



technology to simple models of scalar-driven |l Lj and A-driven |12| inflation. 

An embarrassing postscript to these labors is that the slow roll approx- 
imation purges ^4[^](a;) of its nonlocality. In fact it reduces to —6/R(x), 
which we initially rejected as being dominated by local fluctuations! There 
is actually no obstacle to making use of *4.[g](x) in perturbation theory be- 
cause small fluctuations are, by definition, sub-dominant to the homogeneous 
background. However, one would still prefer an observable which represents 
a more evenly weighted average over the past lightcone. Several alternatives 
are discussed in section 6. We have taken the additional trouble to construct 
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them to reduce exactly to the Hubble constant for homogeneous and isotropic 
geometries, without recourse to the slow roll approximation. 
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